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Abstract 

Security of information has become a major issue during the last decades. 
New algorithms based on chaotic maps were suggested for protection of dif- 
ferent types of multimedia data, especially digital images and videos in this 
period. However, many of them fundamentally were flawed by a lack of ro- 
bustness and security. For getting higher security and higher complexity, in 
the current paper, we introduce a new kind of symmetric key block cipher 
algorithm that is based on tripled chaotic maps. In this algorithm, the uti- 
lization of two coupling parameters, as well as the increased complexity of the 
cryptosystem, make a contribution to the development of cryptosystem with 
higher security. In order to increase the security of the proposed algorithm, 
the size of key space and the computational complexity of the coupling pa- 
rameters should be increased as well. Both the theoretical and experimental 
results state that the proposed algorithm has many capabilities such as ac- 
ceptable speed and complexity in the algorithm due to the existence of two 
coupling parameter and high security. Note that the ciphertext has a flat 
distribution and has the same size as the plaintext. Therefore, it is suitable 
for practical use in secure communications. 

1 Introduction 

Chaos theory is a blanketing theory that covers most aspects of science, hence, it 
shows up everywhere in the world today: mathematics, physics, biology, finance, 
computer and even music. As we know, chaotic systems have many interesting 
features, such as the sensitivity to the initial condition and control parameter, 
ergodicity and mixing property [H [21 El H], which can be connected with some 
cryptographic properties of good ciphers, such as confusion/ diffusion, balance and 
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avalanche property [5l [6]. These properties make the chaotic systems a worthy 
choice for constructing the cryptosystems. Compared with traditional cryptosys- 
tems [7] , the ones based on chaos are more suitable for large-scale data encryption 
such as images, videos or audio data. Furthermore, chaos-based algorithms have 
shown their good performance. Thus, it is a natural idea to use chaos as a new 
source to construct new encryption systems. Cryptography has become one of the 
fields of modern science since the celebrated Shannon's work [8]. In general, the 
algorithms used for cryptographic applications are classified into two, namely, pub- 
lic key (Asymmetric methods) cryptography and secret key (Symmetric methods) 
cryptography. According to the structure of the symmetric algorithm, ciphers can 
be divided into two categories, namely, block ciphers [9] and stream ciphers |10j . 
The chaotic cryptography technique which we are going to use in this paper belongs 
to a symmetric block cipher. There exist two main approaches of designing chaos- 
based cryptosystems: analog and digital. The first one is generally based on the 
concept of chaotic synchronization, initiated by the work of Pecora and Carroll 
and the concerned chaotic systems are implemented in analog form. The second 
one is independent of chaos synchronization and the chaotic systems are completely 
implemented in digital computers [11|T2|. The present paper mainly focuses on the 
digital chaotic ciphers. 

In the digital world, which is currently evolving and changing at such a rapid pace, 
the security of digital information has become increasingly more important. This 
is due to the communications of digital products over network occurring more and 
more frequently. So that in recent years, many different chaotic cryptosystems 
in digital domain have been proposed [131 El 113 [HI 113 HI IIH]- Certainly, a 
main issue in all encryption techniques is their security. However, the recent de- 
velopment of chaotic cryptosystem is rather disappointing. On the other hand, 
various cryptanalysis have exposed some inherent drawbacks of chaotic cryptosys- 
tems [20l HU |22l [23]. Some of these shortcomings are overviewed in |2ll [25l 126] . 
For enhancing the security of discrete chaotic cryptosystem in this paper for the 
first time, we introduce the concept of using hierarchy of tripled chaotic maps with 
an invariant measure in cryptography. Six points make this new cryptosystem 
distinctive and advantageous compared to the other schemes. 

• A very large number of fully developed chaotic maps defined in two intervals 
X G [0, 1] and x G [0, oo). 

• Bifurcation without any period doubling; bifurcation from a stable single 
periodic state to chaotic ones without having usual period doubling or period- 
n-tupling scenario. 

• The existence of the two coupling parameters. One advantage of using two 
coupling parameters is for computational complexity goal. Another advantage 
is the usage of coupling parameters in the structure of the cryptosystem for 
diffusion. These parameters can be used as secret keys as well. 
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• High complexity due to high dimensionahty and chaoticity. 

• Large key space; It is obvious that the attack complexity is determined by 
the size of the key space and the complexity of the verification of each key. 

• The flexibility of attributing different values to the control parameters. 

Besides the above advantages, the encryption speed is also acceptable and the 
ciphertext has the same size as the plaintext. In fact, the results presented in this 
paper totalize our earlier work [27] . 

This paper has been organized as follows: Section 2 describes the construction and 
the main features of the tripled chaotic maps. Current encryption schemes are 
introduced briefly in section 3. The experimental results are shown in section 4. 
Section 5 demonstrates the security of the proposed scheme. Finally, conclusions 
are drawn in Section 6. 



2 Tripled Chaotic Maps 

We first review the one parameter families of trigonometric chaotic maps which 
are used to construct the tripled chaotic maps. One-parameter families of chaotic 
maps of the interval [0, 1] with an invariant measure can be defined as the ratio of 
polynomials of degree N j28j: 

Where F substitute with chebyshev polynomial of type one T]^{x) for 

and chebyshev polynomial of type two Un{x) for a;)p9]. As an example we 

give below some of these maps: 

a^{2x - If ,(2) Aa'^x{l - x) 



,{i) _ a y^x - I) (2) 



4x{l-x) + a^{2x-iy' l + 4(a2_i) 

— x) ' 

The map $2 (x, a) is reduced to logistic one with a = 1. One can show that these 
maps have two interesting properties. It is shown that these maps have interesting 
property, that is, for even values of the (^'^^\a,x){<^^'^\a,x)) maps have only a 
fixed point attractor x = l(x = 0) provided that their parameter belongs to interval 
(iV, 00) ((0, j^)) while, at a > (a > -^) they bifurcate to chaotic regime without 
having any period doubling or period-n-tupling scenario and remain chaotic for 
all a G (0,A^) (a G (-^,cxd)) but for odd values of N, these maps have only fixed 
point attractor at x = for a G {j^,N), again they bifurcate to a chaotic regime 
at a > j^, and remain chaotic for a G (0, -^), finally they bifurcate at a = N to 
have X = 1 as fixed point attractor for all a G (;^,oo). We used their conjugate 
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or isomorphic maps. Conjugacy means that the invertible map h{x) = maps 
/ = [0, 1] into [0, oo) and transform maps $Ar(x, a) into $Ar(x, a) defined as: 

<l>lj'*(a;, a) = — tan^(A^ arctan v^), (2) 

^^^\x , a) = — cot'^ (N arctan —=) , (3) 

Using the hierarchy of famihes of one-parameter chaotic maps (1), we can generate 
new hierarchy of tripled maps with an invariant measure. Hence, by introducing 
a new parameter e as a couphng parameter we form couphng among the above 
mentioned maps, where they can be coupled through (3 and a functions defined as: 

P{x{n)) = {^o + ex{nyf, (4) 

(5) 
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P{x{n + 1)) 



with Bf^{x) and Af^{x) defined as: 

[f] [Vl 
A^ix) = J2Cgx\ Bm{x) = C^k+ix\ (6) 

fc=0 fc=0 

Now, the hierarchy of the tripled chaotic maps can be defined as: 

r Xi{n + 1) = $Arj(xi(n),ai(x2(n),X3(n))), 
^N^,N2,N3 = I X2{n+ 1) = ^N2{x2,a2{x3{n))), (7) 
[ X3(n+ 1) = $Ar3(x3(n),a3), 

Also, their conjugate or isomorphic maps by considering the map h{x) = are 
defined as follows: 

Mn + l) = tan^lA^i arctan Jxi(n)), 



a'({x2{n),X3{n)) 

X2in + 1) = tan2(A^2 arctan y'x2(n)), 

X3{n + 1) = ^ tan^(Ai'3 arctan -I /a;3(n)), 



Ni,N2,N3 



As we know, each dynamical system has specific characteristics. In fact, the dif- 
ferences between discrete dynamical systems arise from these properties. Due to 
this fact, in this paper, tripled chaotic maps are constructed for the first time. So, 
it is necessary to describe them through the paper. Specially, the most important 
characteristics such as invariant measure and Lyapunov characteristic exponents 
have been discussed. 

Furthermore, many properties of the chaotic systems have their corresponding coun- 
terparts in traditional cryptosystems, such as: ergodicity and confusion, sensitivity 
to initial conditions/control parameter and diffusion [30]. For tripled chaotic maps, 
we have tried to describe ergodicity from the invariant measure point of view. 
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2.1 Invariant measure 

Dynamical systems, even apparently simple dynamical systems which are described 
by maps of an interval can display a rich variety of different asymptotic behavior. 
On measure theoretical level these type of behavior are described by SRB [211 [22] 
or invariant measure describing statistically stationary states of the system. The 
probability measure /i on [0, 1] is called an SRB or invariant measure of the map 
y = ^j\f{x,a) given in Eqs. (2)and(3), if it is <I>Ar(a;, a;)-invariant and absolutely 
continuous with respect to Lebesgue measure. For deterministic system such as 
$Ar(x, a)-map, the <l>iv(x, a)-invariance means that, its invariant measure fi{x) ful- 
fills the following formal FP integral equation: 







This is equivalent to: 



^(y)= /.(x)^, (9) 

defining the action of standard FP operator for the map $Ar(x,a;) over a function 
as: 

P'^.fiy)= E (10) 

We see that, the invariant measure fi{x) is actually the eigenstate of the FP operator 
corresponding to largest eigenvalue 1. As it is proved, in our previous work, the 
invariant measure lJ'<i>r^{x,a){x, (3q) has the following form (for more detail see [28]): 

l^ix) = - ^ , (11) 

^y'x(l-x)(/?o + (l-/?o)x) 

with (3q > 0, is the invariant measure of maps Eqs. (2)and(3) provided that, we 
choose the parameter a, in the following form : 

r (iV-i) i 

^^L^=oW^^ (12) 

El. 2 I r~iN Q-k 
k=0^2kP0 

(1 2) 

in $]y' (x, a) maps for odd values of N and 



_ PoJ2 k=0 ^2fc/^0 ^ /iq\ 



2 '/^N O—k 
k=0 ^2fc+ir( 
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in ^^^\x,a.) maps for even values of N. Similarly the probability measure /j, for 
tripled chaotic maps ^Ni,N2,N3 given in Eq. (7) fulfills the following formal FP 
integral equation, where for simplicity we consider their conjugate maps Eq. (8): 

fJ'{yi,y2,y3) ^ j dXij dX2 j dX35{yi- Xi{xi,X2,X3))S{y2- X2{xi,X2)) 

X5{y3 - X3{x3))lJ,{xi,X2, X3), 

which is equivalent to: 

/^(yi,?/2,?/3) = \ J(xi,X2,X2),\ll(xi,X2,X3), (14) 

where J{xi,X2,X2) is the Jacobian of the transformation ^Ni,N2,N3 which is equal 
to: 

d{xi,X2,X3) 
J{Xi,X2,X2) = -TT, T, 

o{yi,y2,y3) 

again defining the action of standard FP operator for the map ^Ni,N2,N3 over a 
function as: 

P^N„N2,N3f^y)^ E f{x)J{Xl,X2,X3), (15) 

^^^N^,N2,N3 

We see that, the invariant measure //(xi, ^2, X3) is also the eigenstate of the FP op- 
erator P^^ ^ ^ corresponding to largest eigenvalue 1 . As it is proved in Appendix 
A the invariant measure /x^^^ (a^i, X2, X3) has the following form: 



1 ^0 1 V/^(^3) 

A* = = X 



^X3{1-X3){(5q + (1 - (5o)x3) ^ ^X2{1-X2){(3{X3) + (1 - (3{x3))x2) 



1 yJP{x2,X3) 
X 1 > (16) 

^ Y'Xi(1 - Xx){l3{x2, X3) + (1 - I3{X2, X3))xi) 

with /? > and given in Eq. (4) and Eqs. (12)and(13). 
2.2 Lyapunov characteristic exponent 

The method of measuring disorder in a dynamical system, is based on the concept 
of Lyapunov Characteristic Exponents (LCE). For one-dimensional maps, LCE 
characterizes the local stretching, determined by the mapping (i0„(a;,a), weighted 
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by the probability of encountering that amount of stretching, that is probabihty 
of a trajectory visiting a particular location x. Here we compute LCE for tripled 
chaotic maps as the characteristic exponent of the rate of average magnification 
of the neighborhood of an arbitrary point tq = (xiq, X2o, 2:30) and it is denoted by 
A(ro), which can be written as |32]: 



, , (9(($ o $ o .... o ($ o $ o .... o $)2. $ o $ o .... o $), 



n— >oo 

k=0 



X(f ) - lim V In I ^3:Af,(xfc,a2) dxN3{xk,a3) , 



where a;fc = $o$o....o<|). A positive LCE implies that two nearby trajectories 
exponentially diverge (at least locally). Negative LCE indicates contraction along 
certain directions, and zero LCE indicates that along the relevant directions there 
is neither expansion nor contraction. The equality of KS-entropy and sum of all 
positive LCE; 

A,>0 

indicates that in chaotic region, this map is ergodic as Birkhoff ergodic theorem 
predicts [33j . 



3 The Encryption and Decryption Procedures 

In this section the framework of our proposed algorithm is described. The proposed 
cryptosystem is a symmetric key block cipher algorithm based on tripled chaotic 
maps. Note that the Lyapunov characteristic exponent is positive, that is, the 
tripled maps are chaotic in nature. According to [M], a possible way to describe 
the key space might be in terms of positive Lyapunov exponents. For this reason we 
have selected all of the control parameters in the chaotic region. A block diagram 
illustrating the complete procedure of the proposed scheme is depicted in Fig. 1. 
First, plaintext is divided into blocks of the same sizes (8-bit)and the blocks are 
transformed into a matrix Mixn, -In matrix M each element of the matrix repre- 
sents one block, and n is the total number of the blocks. The tripled maps are 
iterated using coupling parameters, control parameters and initial conditions. In 
the iteration process first 100 iterations are ignored to avoid transient effects. The 
matrix element Mixj in each round is encrypted as below: 

C^^. =Mixi XOR {xf mod 256), 
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where is a long integer generated using xi, X2 and £3 and some simple mathe- 
matical operators. In each round coupling parameters and x / are regenerated using 
Cixi-i, Xi, X2 and some simple mathematical operators. While the process reaches 
the last element, the elements of matrix Cixn are reversed and Mixn is set equal to 
new matrix Cixn- Then again the process starts from the beginning. In fact, the 
plaintext is encrypted once from the beginning to the end and once from end to the 
beginning. Thus, a very small change in the plaintext will result in a completely 
different ciphcrtcxt. In this cryptosystcm, the process of decryption is completely 
similar to the encryption process. For decryption process, the only difference is 
that the following new relation is used instead of aforesaid relation. 

Mixi=Cixi XOR {xf mod256). 

As an example, the following maps are selected from the chaotic maps(Eq.8) to be 
used in encryption/decryption process. 



2,2,14 



[Xi,X2,Xs) 



^1(^+1) = aj(x.(n).X3(n)) ^^^'(^ ^^^^^^ V ^l(^))' 

X2{n + 1) = tan2(2 arctan ^X2{n)), (18) 

X3{n + 1) = ^ tan^ (14 arctan J X3 (n) ) , 



with 



a2{xs{n)) 



2(3{x3{n)) 



l + p{x3{n))\ (3{xs{n)) 



/3(x3(n + l)) 



P{x3{n)) 



a2 



0L2 



+ ex2,{n))' 



ai{x2{n),X3{n)) 



2p{x2{n),X3{n)) 



p{x2{n + l),X3{n + l))) 



1 + /3(x2(n), 0:3(71)) ^J I3(x2(n),xs(n)) 
P{x2{n),X3{n)) = L hr—^ he'x2{n)f, 



2-ai 



4 Experimental Results 

In this section, we provide some experimental results to illustrate the performance 
of the proposed chaotic cryptosystem. In order to test the efficiency of the proposed 
chaotic cryptographic scheme, we used the scheme in the following files. 

File 1: Text (.txt) file of size 30 720 bytes; 

File 2: Word document (.doc) file of size 210 944 bytes; 

File 3: Executable (.exe) file of size 487 000 bytes; 

File 4: Audio (.mp3) file of size 980 304 bytes; 
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File 5: Image (.bmp) file of size 65 536 bytes; 

File 6: Video clip (.avi) file of size 1 087 430 bytes. 

The encryption and decryption time for all of source files mentioned above are listed 
in Table 1. 

The length of the ciphertext is the same as thought the plain text. This features is 
another most important features of our introduced cryptosystem(See the last col- 
umn of Table 1). Note that, the previous cryptosystems like Baptista-type chaotic 
cryptosystems [35l[36l[37] had almost double-sized ciphertext. 
In order to compare the performance evaluation of the proposed method with the 
previous work [27j, from the security point of view, we focus on the application of 
this method in image encryption. We assume that source image is 65 536 bytes. 
Fig. 2(a) shows the experimental results with Bird BMP image. Fig. 2(b) is its 
encrypted image with the encryption keys mentioned below. According to Eq. (18), 
encryption keys are chosen as follows: 0:3=66, 52=444, xi=445 as initial conditions 
and e=0.2, e'=0.5 as a coupling parameters and finally Q!3=1.5, a2=0.455, ai=0.2 
as control parameters. We have implemented the proposed algorithm using C++ 
programming language and observed the simulation results on a Pentium- IV 2.4 
GHz Celeron D with 256MB RAM and 80 Gb hard-disk capacities. It seems that, 
our encryption time is acceptable compared to that of encryption times mentioned 

5 Security Analysis 

When a new cryptosystem is proposed, it should always be accompanied by some 
security analysis. A good encryption procedure should be robust against all kinds 
of cryptanalytic, statistical and brute-force attacks. Here, some security analysis 
has been performed on the proposed scheme, including some important ones like 
key space analysis, statistical analysis, etc. The security analysis demonstrated the 
high security of the new scheme, as demonstrated in the following. 

5.1 Key space analysis 

A fundamental aspect of every cryptosystem is the key. An algorithm is as secure as 
its key. No matter how strong and well designed the algorithm might be, if the key 
is poorly chosen or the key space is small enough, the cryptosystem will be broken. 
The size of the key space is the number of encrypt ion/ decrypt ion key pairs that are 
available in the cipher system. Apparently, the attack complexity is determined 
by the size of the key space and the complexity of verifying each key. From the 
cryptographical point of view, the size of the key space should not be smaller than 
2^°° to provide a high level of security [HUH]. If the precision 10~^^, the key space 
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size for initial conditions, control parameters and coupling parameters is over than 
2^°°. It seems that the key space is large enough to resist all kinds of brute-force 
attacks. 

In addition, one attempt to describe the dynamics of the current system is by 
providing bifurcation diagram and determining interesting properties of it. In any 
case, the designer of any chaotic cryptosystem should conduct a study of chaotic 
regions of the parameter space from which valid keys, i.e., parameter values leading 
to chaotic behavior, can be chosen. When many parameters are used simultaneously 
as part of the key, the mutual interdependence complicates the task of deciding 
which intervals are suitable. Only keys chosen from the black region of bifurcation 
diagram are suitable enough [l2l HSj. In this cryptosystem, the interval of the 
initial condition is [0, oo). From the bifurcation diagram it is clear that the maps 
are chaotic for any x with respect to control parameters in chaotic region. In Fig. 
3, we have depicted a portion of the bifurcation diagram while N=2. 

As we can see, within the black region, there isn't any periodic windows, so the 
entire black region is suitable for robust keys. 

5.2 Statistical analysis 

In his masterpiece. Shannon [S] said, "It is possible to solve many kinds of ciphers 
by statistical analysis," and therefore, he suggested two methods of diffusion and 
confusion for the purpose of frustrating the powerful statistical analysis. This is 
shown by a test on the histograms of the ciphered images, on the correlations of 
adjacent pixels in the ciphered image and on the distribution of the ciphertext. 

1. Histograms of ciphered images. By taking a (256x256) sized Bird image as a 
plaintext, the histograms of the plaintext and its corresponding ciphertext are as 
follows. It is clear that the encrypted image is confused and cannot be understood. 
Moreover, the histogram of the encrypted image is uniformly distributed, which 
makes statistical attacks very difficult (See Figs. 4(a) and 4(b)). 

2. Correlation of two adjacent pixels. Statistical analysis has been performed on 
the proposed image encryption algorithm by a test on the correlation of adjacent 
pixels in the plain-image and ciphered image. To analyze the correlations of the 
adjacent pixels, we can use Eq. (19) to calculate the correlation coefficients in 
horizontal, vertical and diagonal [T^ 



E{x) is the estimation of mathematical expectations of x, D{x) is the estimation of 
variance of x and cov{x, y) is the estimation of covariance between x and y. where 




(19) 



where 



E(x) = — y^ X. 



D{x) = -J2{x, 
1=1 



E{x)f. 
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X and y are grey-scale values of two adjacent pixels in the image. We randomly 
choose 1000 image pixels in the plain image and the ciphered image respectively 
to calculate the correlation coefficients of the adjacent pixels in horizontal. The 
correlation coefficients of the adjacent pixels in vertical and in diagonal are calcu- 
lated and listed in Table 2 and the distribution is shown in Figs. 5(a) and 5(b). 
It demonstrates that the encryption algorithm has covered up all the characters of 
the plain image and shows good performance of balanced 0-1 ratio. 
In addition to analysis mentioned above, we have also analyzed the distribution of 
the ciphertext. Therefore, we have recorded the number of occurrences of each ci- 
phertext block for one of the six source files. A typical distribution of the ciphertext 
is shown in Figs. 6(a) and 6(b), which shows that the distribution is very flat due 
to the masking operation. Totally, statistical analysis has been performed on the 
proposed encryption algorithm, demonstrating its superior confusion and diffusion 
properties which strongly resist statistical attacks. 

5.3 Information entropy 

Information theory is a mathematical theory of data communication and storage 
founded in 1949 by Claude E. Shannon. To calculate the entropy H(s) of a source 
s, we have: 



where P{si) represents the probability of symbol Sj. Actually, given that a real 
information source seldom transmits random messages, in general, the entropy value 
of the source is smaller than the ideal one. However, when these messages are 
encrypted, their entropy should ideally be 8. If the output of such a cipher emits 
symbols with an entropy of less than 8, then there exists a predictability which 
threatens its security. We have calculated the information entropy for encrypted 
image Fig. 2(b): 



The obtained value is very close to the theoretical value 8. Apparently, comparing 
it with the other existing algorithms, such as [15], the proposed algorithm is much 
more closer to the ideal situation. This means that information leakage in the 
encryption process is negligible, and so the encryption system is secure upon the 
entropy attack. 

5.4 Differential attack 

As we know, if one minor change in the plain-image can cause a significant change 
in the ciphered- image, with respect to both diffusion and confusion, then this "dif- 




(20) 




255 
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ferential attack" may become inefficient. To test the influence of one-pixel cliange 
on tlie whole image, encrypted by the proposed chaos-based algorithm, two com- 
mon measures may be used: Number of Pixels Change Rate (NPCR) and Unified 
Average Changing Intensity (UACI) [IHl HE] . We take two encrypted images, Ci 
and C2, whose corresponding original images have only one-pixel difference. We 
label the grey scale values of the pixels at grid (i,j) of Ci and C2 by Ci{i,j) and 
C2{i,j), and Ci and C2 have the same size. Then, D{i,j) is determined by Ci{i,j) 
and C2(i,i), that is, if Ci{i,j) = C2{i,j), then, D{i,j) = 1; otherwise, j) = 0. 
NPCR and UACI are defined by the following formulas: 



Where, W and H are the width and length of the image. We obtained NPCR=0.34 
% and UACI=0.33 %. With regard to obtained results, apparently, the proposed 
algorithm has a good ability to resist differential attack. 

6 Conclusion 

In this paper, after introducing hierarchy of tripled chaotic maps with invariant 
measure, we investigated their potential for exploitation as a cryptosystem. These 
maps have interesting features such as invariant measure, ergodicity and variable 
control parameters. In addition to some features aforesaid, the most important 
advantage of these maps is the existence of the two coupling parameters. It seems 
that the excellent efficiency of the new cryptosystem is derived from this property. 
In this scheme, its structural parameters and initial values can all be used as en- 
cryption key in chaotic cryptosystem. Therefore, the size of key space is very large. 
Statistical analysis performed on the proposed encryption algorithm, demonstrates 
its superior confusion and diffusion properties which strongly resist statistical at- 
tacks. Experimental results illustrate that the distribution of the ciphertext is very 
fiat and the entropy, is almost equal to the ideal value. Moreover, the encryption 
time is acceptable while the size of ciphertext is the same as that of the plaintext. 
The present paper's goal is to construct a practical and fully secure cryptosystem. 
Our experiments indicate that this goal is almost achieved. Based on all analy- 
sis and experimental results, the conclusion is that, from a cryptographical point 
of view, the proposed scheme is a good candidates for practical applications in 
information security fields. 



NPCR 




X 100% 




X 100% 
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Appendix A: Derivation of the invariant measure In order to prove that 
measure Eq. (16) satisfies Eq. (14), it is rather convenient to consider the conjugate 
map. By inverting the conjugate map Eq. (8) we get: 



' Xkz = $^(a;3,Q;3)), 

< Xk2,k2. = ^^A^2-, Ol2{Xkz)), 

Now, considering the following anzatz for invariant measure 

1 1{X2:X3.) 



(-1) 



IJ,{Xi,X2,X3 



where 



TT {1 + P{x2,X3)xi)^/x]^' 

1{X2,X3) 



(•2) 



^J(5{X2,X3) 



dX2dX3 



The probability measure for coupled chaotic maps ^Ni,N2,N3 given in Eq. (7) 
fulfills the following formal FP integral equation (Eq. 14): 



fx{Xi,X2,X3) = Yl 



dXkl,k2,k3 dXk2,k3 dXkS , . ^ 

l^[Xkl,k2,k3, Xk2,k3, Xk3) 



dxi dx2 dx3 

by taking the derivative of last term of Eq. (A.l) with respect to Xi, we have: 

dXkl,k2,k3 _ Q;i(Xfc2,fc3, Xks) i/a^fcl,fc2,fc3(l + Xkl,k2,k3) 

dxi ~ Ni 
from Eq. (A.2) it follows that: 

7(^:^2, a;3) 



(.3) 



E 



Xi{l + al{Xk2,k3,Xk3)xi) 
(1 + Xkl,k2,k3) 



;-4) 



(1 + I3{X2, X3)xi)^/Xl klMM ^1 + al{Xk2,k3, Xk3)xi) 

l{Xk2,k3,Xk3) 



X 



1 + (3{Xk2,k3, Xk3)Xkl,k2,k3 



(.5) 



we readily see that 

E 



«1 (1 + XklM2M3) _ (^iJ^nAP Kxk2,k3,Xk3)) 
k3 ^1 (1 + I3i.^kl,k2,k3)) 5jVi(/3-n^fe2,fc3, Xk3)) 

{1 + al{Xk2,k3^Xk3))xi) 



X 



1 + l3{Xk2,k3, Xk3){^^^jP^fxi 



:-6) 



Eq. (A. 3) reduce to: 



1{X2,X3) 



[I + (3{x2,X3)xi] 



E 



dXk2,k3 dXk3 

dx2 dx3 
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which is possible if P{x2, X3) and P{xk2,k3, Xks) are related as: 

\ BnAP {Xk2,k3,Xk3)) J 

therefore, the relation Eq. (A. 7) reduces to 

7(2^2,2^3) ^ dXk2,k3dXk3 l{Xk2,k3,Xk3) , . . „s 



^J/3{X2,X3) k2,k3 5a;3 ^P{Xk2,k3,Xk3)Xk3 

By considering 

7(2^2, 2^3) 1 7(3^3) 



^JP{X2,X3) TT ^{1 + P{X3)X2) 

where 

^^"^^ =ciX3 = 1 



which is possible if /^(xa) and P{xk3) are related as: 

/ a2ix3)AN^{p-^{xk3)) 



or Eq. (A. 9) can be written as: 



1{X3) Y^dxk3 ^{xks) 1 (3 



.10) 



/3(Xk3) = /3(X3) (.11) 



/P{X3) ^ 9X3 ^l3{Xk3) ^ V^(l + /^^3) 

hence it is proportional to its invariant measure. 
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Table 1: Performance of the proposed chaotic cryptographic on the six files, en- 
cryption / decryption t ime 







Plaintext file size 


Encryption time (in seconds) 


Ciphertext file size 






(bytes) 


min.-max. (mean) 


(bytes) 


File 1 


(.txt) 


30 720 


0.14-0.15 (0.145) 


30 720 


File 2 


(.doc) 


210 944 


1.005-1.007 (1.006) 


210 944 


File 3 


(.exe) 


487 000 


2.5-2.56 (2.53) 


487 000 


File 4 ( 


.mp3) 


980 304 


0.44-0.47 (0.455) 


980 304 


File 5 ( 


.bmp) 


65 536 


0.33-0.36 (0.345) 


65 536 


File 6 


(.avi) 


1 087 430 


5.61-5.67 (5.64) 


1 087 430 
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Table 2: Correlation coefficients of two adjacent pixels in two images 





Plain image 


Ciphered image 


Horizontal 


0.9902 


0.000029549 


Vertical 


0.9815 


0.00045066 


Diagonal 


0.9721 


0.000013 
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Figures Caption: 
Fig. 1: Block Diagram. 
Fig. 2: (a) Plain image, (b) Ciphered image. 
Fig. 3: Bifurcation diagram of ^^^\x,a) while N=2. 
Fig. 4: (a) Histogram of plain image, (b) Histogram of ciphered image. 
Fig. 5: (a) Correlation analysis of plain image, (b) Correlation analysis of 
ciphered image. 

Fig. 6: (a) Distribution of the plaintext, (b) Distribution of the ciphertext. 
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